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}Q ! Abstract 
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A number of 2d and 3d four-fermion models which are renormalizable — in the 1/N 
expansion — in a maximally symmetric constant curvature space, are investigated. To this 
purpose, a powerful method for the exact study of spinor heat kernels and propagators on 



maximally symmetric spaces is reviewed. The renormalized effective potential is found for 
any value of the curvature and its asymptotic expansion is given explicitly, both for small and 
for strong curvature. The influence of gravity on the dynamical symmetry breaking pattern 
of some U(2) flavor- like and discrete symmetries is described in detail. The phase diagram 
in S 2 is constructed and it is shown that, for any value of the coupling constant, a curvature 
exists above which chiral symmetry is restored. For the case of H 2 , chiral symmetry is always 
broken. In three dimensions, in the case of positive curvature, S 3 , it is seen that curvature 
can induce a second-order phase transition. For if 3 the configuration given by the auxiliary 
fields equated to zero is not a solution of the gap equation. 
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1 Introduction 



Four-fermion models JT], |2| — usually considered in the 1/N expansion — are interesting due 
to the fact that they provide the opportunity to carry out an explicit, analytical study 
of composite bound states and dynamical chiral symmetry breaking. At the same time, 
these theories — and specially their renormalizable 2d and 3d variants — exhibit spe- 
cific properties which are similar to the basic behaviors of some realistic models of particle 
physics. Moreover, this class of theories can be used for the description of the standard 
model (SM) itself, or of some particle physics phenomena in the SM (see For ex- 

ample, the dynamical symmetry breking pattern of Nambu-Jona-Lasinio (NJL) models for 
the electroweak interaction, with the top quark as an order parameter, has been discussed 
in SB- 

Having in mind the applications of four-fermion models to the early universe and, in 
particular, the chiral symmetry phase transitions that take place under the action of the 
external gravitational field, there has been recently some activity in the study of 2d ||, 
3d |J and 4d ||10||-||14| four-fermion models in curved spacetime (for a general introduction 
to quantum field theory in a curved spacetime, see |15|)- The block-spin renormalization 
group (RG) approach and the similarities of the model with the Higgs- Yukawa one have been 
considered in [TB| and WM, respectively. 

The effective potential of composite fermions in curved spacetime has been calculated in 



different dimensions fl8|-||ll|. Dynamical chiral symmetry breaking, fermionic mass genera- 
tion and curvature-induced phase transitions have been investigated in full detail. However, 
in most of these cases only the linear curvature terms of the effective potential have been 
taken into account |9[|-[1l^|. But it turns out in practice that it is often necessary to con- 
sider precisely the strong curvature effects to dynamical symmetry breaking. In fact we will 
see that going beyond the linear-curvature approximation can lead to qualitatively different 
results. 

In this paper we will investigate some 2d and 3d four-fermion models which are renormal- 
izable — in the 1/N expansion — in a maximally symmetric constant-curvature space (either 
of positive or of negative curvature). The renormalized effective potential will be found for 
any value of the curvature and the possibility of dynamical symmetry breaking in a curved 
spacetime will be carefully explored. Furthermore, the phase structure of the theory will be 
described in detail. 

The paper is organized as follows. In the next section we calculate the effective potential 
of composite fermions in the Gross-Neveu model, in the spaces S 2 and H 2 . The phase 
diagram in S 2 is constructed and it is shown that for any value of the coupling constant there 
exists a curvature above which chiral symmetry is restored. For the case of H 2 , we show 
that chiral symmetry is always broken. The asymptotic expansions of the effective potential 
are given explicitly, both for small and for strong curvature. The three-dimensional case 
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is studied in Sect. 3. We consider two different four-fermion models: one which exhibits 
a continuous U(2) symmetry and another where we concentrate ourselves on two discrete 
symmetries which happen never to be simultaneously broken (see In Sect. 4 we study 

explicitly the dynamical P and Z 2 symmetry breaking pattern in H 3 and S" 3 . Finally, Sect. 
5 is devoted to conclusions and some technical points of the procedure are summarized in 
two appendices. 



2 The Gross-Neveu model in a space of constant cur- 
vature 



2.1 Case of the 2d De Sitter space S 2 

In this subsection we will undertake the discussion of the Gross-Neveu model in De Sitter 
space. This model, although rather simple in its conception, displays a quite rich structure, 
similar to that of realistic four- dimensional theories — as renormalizability, asymptotic free- 
dom ||2|, [I7j and dynamical chiral symmetry breaking. Some discussions of chiral symmetry 
restoration in the Gross-Neveu model for different external conditions (such as an electro- 
magnetic field, non-zero temperature or a change of the fermionic number density) have 
appeared in the past |TB|, |2T], [3D] (the influence of kink-antikink configurations on the phase 



transitions is described in 22 



The study of the Gross-Neveu model in an external gravitational field has been performed 
in Ref. || using the Schwinger method [E4l] (for other analysis of two-dimensional models 



in curved space, see f25|). Unfortunately, the generalization of the Schwinger procedure to 
curved spacetime is not free from ambiguities and this is why the result of Ref. || includes 
some mistake. 

In the calculation of the effective potential in the Gross-Neveu model on De Sitter space 
we will use a rigorous mathematical treatment of the fermionic propagator in (constant 
curvature) spacetime, that has been developed in Ref. [p6| . We shall start from the action 



S 



d 2 x y/—g 



(i) 



where N is the number of fermions, A the coupling constant, 7 M (x) = 7 a e^(x), with 7° the 
ordinary Dirac matrix in flat space, and V M is the covariant derivative. By introducing the 
auxiliary field a, it is convenient to rewrite ([[]) as 



S 



d x 



(2) 



with a = —j^ipip. Furthermore, in order to apply the results of Ref. |[26| , it is convenient to 
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use Euclidean notations. Then fl2T) is written as 



S 



ip^(x)W^ + ^a 2 + atpip 



(3) 



Assuming that we work in De Sitter space and using the standard 1/N expansion, we 

tern 



get the effective potential in terms of the a field as follows || [10], |TT 

■2 



a 
2A 



Tr / D(x,x,s)ds, 
Jo 



where the propagator D is defined by 

(V + s)D(x, y, s) = -S 2 (x, y). 

The curvature of both S and H d can be written in the form 

d(d-l)k 
R — 5 , 

a 2 

with k = 1 for S d and k = — 1 for H d ; a stands for the radius of the manifold. 



(4) 



(5) 



(6) 



We consider first k — 1. Following we begin the calculation by obtaining the 'square- 
ing' Green's function: 

(v 2 - s 2 )G(x,y) = -5 2 (x, 2/ ). 
Then D = (V — s) G(x,y). We resort to the Ansatz p6| 



(7) 



G(x,y) = u{x,y)g(p), (8) 

where u(x,y) yields a unit matrix in spinor indices when y — > x, and p is the distance 
between x and y along the geodesic that goes through these two points. Introducing the 
notations 

= Z g(6) = cos 9 - h(6), Z = cos 29 -, (9) 
this leads to the following equation for h 

d 2 



Z(l-Z) 



dZ 2 



(2-3Z)A-l-sV 



h{9) = 0. 



(10) 



One can show [2!J that the solution of (^) is given by a linear combination of the hypergeo- 
metric functions |27| . As boundary conditions for (^) it is convenient to choose the following. 
First, one selects the singularity that appears for 9 — > 0. Second, one demands [26] that the 
singular part of this limit has the same form as in flat space 



g (6)~- — ln0, 6^0. 

Z7T 



(11) 



Using properties of the hypergeometric functions and the boundary conditions (|TTD, the 
function D for coinciding arguments is found to be 



D(x, x, s) 



-s lim q(9) = — 



9 2 

ip(l + isa) + ^(1 - isa) + 27 + In — 



(12) 
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where (|TT|) has been used explicitly, 7 is the Euler constant and ip the digamma function, and 
the two arguments of D are supposed to be separated by a small geodesic distance p = 9a. 
Differentiating (H) with respect to a and using (JT3), one gets 



^ + TYD(«,,, ( ,) = ^1 + A 



27 + ^ ( 1 + i|o-|^/ — 



(13) 



where R is the curvature. 

As renormalization conditions we choose the following 



V"{a) 

By selecting the counterterms of the form 

5V = 



i?=0,cr=/x \ ' 



1 2 /ipe 7+1 

— cr In 

2tt 2 



(14) 



(15) 



adding them to fli~3|) and using the asymptotics of i/j(x), we find the following value for the 
derivative of the renormalized effective potential 



V'(a) 



a 
A 



1 



2^ 



2 - In 



(16) 



As we can see, letting aside the different notation employed for a, the terms involving ijj(x) 
differ from the corresponding terms reported in |J. 

Starting now from expression (|16|) , different physical questions can be studied. In par- 
ticular, the possibility to construct a corresponding phase diagram appears. To this end, let 
us calculate the second derivative 



V"{a) 



1 1 
A + 2~7r 

i\a\ 



' ' 1 "" '^^r \ + ^ ( 1_i ' (T 'vl' 



2 - In 



2tt V R 



C |2,l + iH^/|j -c [2,1-W| 



(17) 



We can now study the behavior of the renormalized V near a = 0. We have always 



no) = 0, v"(o) = j 



(18) 



(19) 



With the notation 

Rq = 2/iV (7+1) , A = 27r, 

we obtain that the point a = is a minimum for R > i?o exp(— Ao/A) and a maximum 
for R < Roex.p(—\ /\). That is, for any value of A, there exists a value of the curvature 
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above which chiral symmetry is restored. The connection of the critical curvature with the 
coupling constant is 

R cr = Roe- Xo/x , (20) 

a different expression for R cr will be given below, in which the independence of R cr from the 
renormalization scale /i is made apparent. For R < R cr the chiral symmetry is broken and 
a dynamical fermion mass is generated. 

As next step one can investigate different limits of the expression (|16"D for the renormalized 
potential. For R — > 0, one obtains from (16) 

- R ' 



Via) 



a 
A 



1 



A 
2^ 
A 



a 



In — - 2 + 



12a 2 



1 + — In — - 3 



2tt 



R 



48tt R 



21 



(22) 



Notice that in (|2 
of V as R -> 



the constant of integration has been chosen having in mind the finiteness 
At the same time, at R — ► oo one gets 

27 
A 



V'fcr 



(7 

A 



1 + 



-27 



a- 



2tt 



a 2 















(23) 
(24) 



The last expressions show the behavior of the effective potential at strong curvature. Analysing 
(|2T|)-(|24D one can see that for small R chiral symmetry is broken, as it happens in flat space- 
time. However, in the limit of strong curvature chiral symmetry is restored. Thus, the study 
of the asymptotics is a further check of our general analysis (phase diagrams). 
If one defines M by 

V"{a)\ R=0>a=M =0, 

then one can express the derivative of the effective potential in terms of this parameter, as 

R 



a 



In 



^1+14 



~R 



+ V 1 



12a 2 
~R 



In 



2M 2 



We may now rephrase the criterion of symmetry restoration in terms of M by saying that 
the symmetry is restored when R > R cr = 2M 2 exp (27). The shape of the effective potential 
for different values of the quotient is shown in Fig. 1. The character of the transition is 
continuous, as illustrated by Fig. 2. 



2.2 Case of the hyperbolic space H 2 

In the hyperbolic space H 2 (negative curvature) the analysis can be carried out in a very 
similar way. After introducing the notations 



9 



P 
a 



6 9 
g{9) = cosh - h(9), Z = cosh 2 — . 



(25) 
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Eq. (0) acquires the form 



Z(l - Z)^L + (2 _ 3 Z) A _ , + sV 



As a result, using a similar procedure as in Sect. 2.1, we get 



h{9) = 0. 



(26) 



x, s) 



Air 



27 + + sa) + if)(sa) + In 



0' 



(27) 



The non-renormalized expression for V'(a) is 



^+*( 1 +H/if)+*(H/ii)+ h?; r 



(28) 



Making the same renormalization as in Sect. 2.1, we obtain the renormalized effective 
potential 



V'{a) 
V"(a) 



a f A 
a{ 1 + 27 

1 I, A 

A I 27T 



2 H 1 + H \/M)- 2 - ln lT 



8tt 2 



(29) 



(30) 

A careful study of V"(0) shows that due to the presence of the last term in (|2"9"|), a = is 
never stationary for any value of A and finite R. Owing to the fact that V"(0) < 0, chiral 
symmetry is always broken in H 2 . 

In the small curvature limit (\R\ — ► 0) 



V(a) 



2A 



\R\ , a 2 
m 



48vr \R\ 



(31) 



what coincides with fl22|) , taking into account the change of sign for the curvature. In the 
strong curvature limit (\R\ — ■> 00) 



cr 



(32) 



The analysis of Eqs. (31]) and (|3^) shows that the general conclusion about the chiral 
symmetry breaking at any finite R in H 2 is correct. In a similar way one can study the 
influence of curvature in the massive Gross-Neveu model (for a recent discussion of such 
model at non-zero temperature, see fl2T 
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3 Dynamical U(2) flavor symmetry breaking in H 3 and 

s 3 

This section is devoted to the description of a three-dimensional four-fermion model which 
has a continuous flavor-like symmetry and how its breaking is affected by a gravitational 
background. 

We consider the following model on a Riemannian manifold 



2N 



(33) 



with 



r 



^ 

-a* 



, a* = 1,2,3, T 



1 2 

-i 2 o 



-AV7 3 , 



1 2 

i 2 



(34) 



so we take a reducible, four- dimensional Dirac algebra. 

The transformation of the bilinear terms ipifj, ifjir 1 ^ and ipi^ip under 

with T = 1, Ti = 7 5 , T 2 = r 1 and T 3 = frV, is given by 



(35) 







/ 


-0 1 -9 2 \ 


( i>i> \ 




= 2 




-# 3 




v 5(^7^) y 




V # 2 


3 j 





(36) 



It is convenient to express this theory in terms of auxiliary fields, namely 

C E = i) (0 + a + iT X p + i 7 5 vr) V + ^(^ 2 + P 2 + 7T 2 ), (37) 

what yields 

f N 

T N _ 00 [a,p,iT} = J dx^ — ( ( T 2 + p 2 + 7r 2 )-A^lndet(^ + ( j + ir 1 p + «7 5 7r)- (38) 



We consider the constant configurations of a, p, n. The rotational symetries impose that the 
effective potential only depends on a 2 + p 2 + it 2 , what allows us to set p = it = 0. The 
regularized expression of the effective potential turns out to be: 



r r - N n N r 

J dx^gV[<j]=J dxy/g— a 2 + -Tr ^ 



<-°° dt 

/A 2 t 



-ta 2 V 2 t 



(39) 



The reader may verify that the use of point-splitting regularization leads to the same results 
that we get below by regulating the lower limit in the proper time integral. On the right 
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hand side we observe the appearance of the coincidence limit of the heat kernel. Resorting 
to the results of appendix B, the outcome is (for H 3 ) 

Via] = -—cr 2 H ttt / -tat ( H )e ta e y \ (40) 

1 J 2X B 2tt 3 / 2 A/A2 t 5 / 2 V 2tJ K ' 

The dependence on the radius of the manifold is not shown, but it can very easily be recovered 
from dimensional analysis. We obtain 















+ 2aVr (~, 


-11 




a 2 ; 




A 2 J 



N 2X B (4vr) 3 / 2 a 2 

Imposing the renormalization condition 

1 d 2 _ . , , I 



(41) 



V[a = 0,R= -Q/a 2 = 0] = -, (42) 



one gets 



and in the limit A — > oo: 



Nda 2 1 ' ' ! A 

1 1 A 



\ B A 7T 3 / 2 ' 



(43) 



N V W = TX + i{-3-^)- (44) 

From this one can easily see that the symmetry is always broken in H 3 . In fact we can also 
notice the curious feature that the origin is not a solution of the gap equation. This appears 
to liken the situation of four-fermion models in three dimensions under the influence of a 
magnetic field (see p8| ). 

As for S 3 , we have (see appendix B). 

a 2 1 r 2 1 

V{°) = r-— rTr / d s . (45) 

v ; 2\ B 2 Vol Jo s- p y ' 



From this formula one derives the regularized effective potential using the result from [26 
that on S N the solution of 

(W 2 -™ 2 )G(y) = -5 N (y) 

is given by G(y) = U(y)gi>f(a), where U(y) is a parallel transport matrix that is the identity 
at the coincidence limit, and 

r(f + im)r(f -im) (e\^(N . n . n , 2 (e\\ 

g N (0) = — ^ L - T rr — 7-^ cos - F — + im, im, — + 1, cos 2 - , 

KJ ( 47r )^ 2 r (f + 1) V 2 / V 2 2 2 V 2 // 

where 9 = E , (p is the geodesic distance, and a is the radius of the manifold). With this 
in mind and using known properties of the hypergeometric functions one may arrive at the 
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following expression (in which we only keep terms divergent in p and terms independent 
from it) 

tw/ x ° d Tr 1 4 Tr 1 4 / 1 9 \ , . . , , . 

v ' \ B An p An \Aa 2 J 

Using the renormalization given by 

i^ = 0,fl = 6/«' = 0] = i, 

we get 

V\a) = ^ + - ( -L + a 2 ) tanh (ixoa) . (47) 

This appears to be in agreement with the correspondent effective potential in the Gross- 
Neveu model (in the simplest version without continuous U (2) symmetry and discrete sym- 
metries) on De Sitter space using, different from our approach, dimensional regularization 
(see We can now compare this result with the one that was found in Ref. f|, which 

was a study of three-dimensional theories in the small curvature limit. One can easily check 
that expression (|46f) is indeed compatible with those in ||. What is most surprising in the 
weak curvature limit is that, in view of the results, one might be tempted to conclude that 
the origin is not a solution to the gap equation (as it happens in H 3 ), and, furthermore, that 
there may be no solution at all. However, looking at the exact result fl46|) or (0), this is 
seen to be an artifact of the approximation. The effective potential given in expression (|47D 
may give rise to a second order phase transition, as we illustrate in Fig 3. 

4 Dynamical P and Z2 symmetry breaking in H 3 and 

s 3 

In this section we analyse a model which displays two discrete symmetries. First of all, we 
present the model and later we describe its symmetries in some detail. To finish, we will 
describe the influence of gravity on the breaking of the symmetries. Using the representation 
for the 7 M (which has no 75) 

it is also possible in this case to define a parity operation which admits the presence of a 
mass term, i.e. 

C E = $ jty + mijjT 3 i/j } (49) 

with 

P[ip(x 1 ,x 2 ,x 3 )] = 7 2 tV(xi, -x 2 ,x 3 ), P[i)(xi,X2,x 3 )] = -i>(xx, -x 2 ,x 3 )j 1 t 1 ]. (50) 
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However, ifjip transforms as a pseudoscalar. There is another discrete symmetry of the kinetic 
term, given by 

Z 2 [^]=rV, Z 2 [$\=$t 1 . (51) 
Under this operation, the roles of the mass terms are reversed, e.g. 

Z 2 H>T 3 ij] = -^jr 3 ij, Z 2 [^\ = (52) 

With this in mind, it is immediate that 

C E = $W- ^f(^) 2 - H (^) 2 (53) 

is invariant under both P and Z 2 . But, as iJjt 3 i/j is not invariant under Z 2 , it can be taken as 
an order parameter for the Z 2 symmetry breaking. Likewise, < Tptp > is an order parameter 
for the P symmetry breaking. 
In terms of auxiliary fields 

C E = + + X r 3 )^ + + ^X 2 - (54) 

Proceeding along the same lines we have trodden before, we will impose equivalent renor- 
malization conditions 

<P ; V[<P = 0,R = 0] = \, ^-V[ X = 0,R = 0] = -. (55) 



d<p 2 A dx 2 k 

In the case of hyperbolic space if 3 we get that the effective potential in the large- N limit is 
given by 

1 . o\ + a 2 _ fl IN a + a_ fl IN 

<± \ A H/ 
\ f „1 1 \"1 

(56) 



N 8 VA k 4 V A k 



1 

+ 2i 



with cr ± = cj) ± x- We obtain the two cases 



\-->0, 0+ = -<7-, (0 = O,x^O), (57) 

A K 

for the Z 2 symmetry breaking, and 

\--<0, <J + = a_, (0^O,x = O), (58) 

A K 

for the P symmetry breaking, respectively. This is illustrated in Figs. 4 and 5, where we 
only consider positive values of 0, as the symmetry of the model allows us to reproduce the 
result for the region of negative immediately. It is also worth noting that either P or Z 2 
is broken, but that it is impossible to have both symmetries broken. Fig. 4 exemplifies the 
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first situation (2.68): one sees that the minimum lies at = 0, x 7^ an d that Z2 is broken. 
In Fig. 5 we are in the second situation: the minimum lies at 7^ 0, x = and P is broken. 

As for the S 3 case, the reader may verify — along the same lines as above — that one 
encounters again two different situations: either both P and Z 2 are unbroken or just one of 
them is broken. This is obtained without difficulty by repeating the same analysis, and we 
feel that further details are not necessary. 



5 Conclusions 

In this paper we have reviewed a powerful method for the exact study of spinor heat kernels 
and propagators on maximally symmetric spaces. We have used it with success in a number 
of different four-fermion models. The renormalized effective potential has been found, in each 
case, for any value of the curvature, and its asymptotic expansion has been given explicitly, 
both for small and for strong curvature. 

We have described in great detail the influence of gravity on the dynamical symmetry 
breaking pattern of some U(2) flavor-like and discrete symmetries. In particular, we have 
seen explicitly that the effect of a negative curvature is similar to that of the presence of a 
magnetic field. 

In the case of positive curvature, S 3 , we have checked that the scenario given in the 
framework of a small curvature expansion is dramatically changed when we treat gravity 
exactly. In particular we find that the character of the phase transitions induced by gravity 
is continuous and that the origin is always a solution of the gap equation. A point to be 
duely remarked is the fact that the techniques applied here to two- and three-dimensional 
models work equally well in four- and higher-dimensional ones, as well. Of course, these 
models are not renormalizable in the standard way (see however, refs. |!2| where mean field 



renormalization of four-fermion models has been discussed), but one can still apply the above 
method in order to obtain the cut-off dependent effective potential, what can certainly be 
useful for cosmological applications. 

Another quite interesting topic is to study a formulation overlapping technique between 
the Schwinger- Dyson equations in a constant-curvature space and the effective potential 
approach. We hope to address some of these questions in the near future. 
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A Appendix: Generalities about spinors in three di- 
mensions 



Some general facts about spinors in three dimensions that may be interesting to recall for the 
benefit of the reader are the following. The Dirac algebra {7^, 7^} = 2g fiv has two irreducible 
representations, namely 

1. 7 = a 3 , 7 1 = ia 1 , 7 2 = ia 2 
and 

2- 7; = -7,. 

On the other hand, there is no '7 s ' matrix in order 2 which anticommutes with all of the 
7 M s. A mass term of the form rrnpip in the Lagrangian explicitly violates parity, defined by 

P : ii(x°, x 1 , x 2 ) — ► trty(a: , -x 1 , x 2 ). (59) 

However, if one uses a reducible representation of the Dirac matrices, such as 

'-(0 (60) 

then the mass term m'ipip does preserve parity, as defined by 

P : ip(x°, x\x 2 ) — ► rSV^ , ~x\x 2 ), (61) 

being 

j / 1 2 \ 2 ( \ 3 / 1 2 






1 2 \ 


-( 


V i 2 


' 





The kinetic term is invariant under the U(2) transformations with generators T = 
1 4 , Ti = 7 5 , T 2 = r 1 , and T 3 = ir^ 5 , where 

7 5 = AVV = if ^ (63) 
The mass term breaks this symmetry down to U(l) x U(l) with the generators To, T3. 



B Appendix: Spinor heat-kernel in maximally sym- 
metric spaces (the intertwining method) 

In this appendix we review a powerful method to derive closed expressions for the spinor 
heat-kernel and propagators in maximally symmetric spaces (for more details see p6|). 
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The spinor heat kernel satisfies (— X 72 ) 2/o; t) = and the initial condition lim^o 
K(y,y ;t) = 1 SN(y,yo)- Substituting the Ansatz K(y,t) = U(y)f(cr,t) in the heat-kernel 
equation, one gets 



u % + uvaVaf + 2 ^(v a c/)^ + (V a V a f/)/ - |c// = 0, 



(64) 



where n a = V a o\ Here the term linear in df /da and V ' a U cancels out, provided that U 
satisfies the parallel transport equation 



n a V a U = 0, U(y ) = l. 
The Laplacian acting on / can be replaced by its radial part given by 



V a V a f = a N f = {dl + (N — l)Bd a )f, B 



cot 



S N , 

±cothf^, H N , 



and the Laplacian acting on U is 

V a V a U = ~A 2 (N - 1)U, 

The equation for the scalar / becomes 

{-d t + L N )f = 0, L N 



itanh(^), H», 



□at- | -\{N-l)A\ 



(65) 



(66) 



(67) 



(68) 



The idea is now to relate the solutions of this equation for different iVs. To this end one 
looks for an operator D such that L N D = DL N _ 2 - A simple Ansatz leads to 

«' ^ ' (69) 



D 



JL C o S I 9 
2tt LUb 2 d(cos 0) 

— ± cosh | 



cos 



i) ( cosh f) , 



X 



H 



N 



The odd-dimensional case is elementary, taking into account that L^D^ N = D^ N l ^ 2 L\ 
and the known solution of (—d t + Li)K\ = 0. This yields 



K N (y,t) 



x 



U (y) cosh - 



d 



(JV-l)/2 



2ir <9(cosh:r) / 

oo 

K N (y,t) = U(y) ]T (-l) n f dp (0 + 2nn,t), 



( cosh — 
V 2 



-1 e -^ 2 /(4t) 



H 



N 



(70) 



/c«p(^ + 27i~n, t) = cosh 



( 1 



9 



2 ^Tr^coshfl)) 



(JV-l)/2 



cosh 



-i e -e 2 /(4t) 



In the particular case N = 3 these expressions read — in the coincidence limit — as follows 



lim K(y,t) = t d — 
y^yo (Anty/ 2 



1 + 



t 

2^2 



H 3 



1 



2^2 



+ £(-!)» exp 



n=l 



7r 2 a 2 n 2> 
t 



lim = l d , 1- 

y^2/o (47rt) 3 / 2 

where d is the dimension of the Dirac algebra, in our case d — 4. 



2 - t - 4vr z n 



2„2' 
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Figure captions 

Fig. 1. In this figure, v = Va 2 , s = era and r = j^. One sees that for r > rc = jfi the 
symmetry is restored in S* 2 . 

Fig. 2. In this figure, r = and min = aa min , where <7 m j„ is the value of <r at which V 
attains its minimum. The continuous character of the transition which takes place in S 2 is 
clearly exhibited. 

Fig. 3. Here, v = Va 3 , s = aa and g = —. It is clearly seen that there may be a second 
order phase transition in S* 3 and that it takes place when g — 4. 

Fig. 4. This figure is a plot of Va 3 in i^ 3 in terms of s = 0a and t = We have taken 
f = 1 and * = 0.25, so j - 1 > and Z 2 is broken. 

Fig. 5. This is also a plot of Va 3 in if 3 , with s and t defined as in Fig. 4. Here we have 
y — - < and P is broken. 
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Figure 4 




Figure 5 




